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SENSITIVITY ANALYSIS OF TORSIONAL" VIBRATION 
CHARACTERISTICS OF HELICOPTER ROTOR BLADES 


PART I - STRUCTURAL DYNAMICS ANALYSIS 

By Theodore Bratanow* and Akin Ecer** 

University of Wisconsin-Milwaukee 

SUMMARY 

A theoretical investigation of structural vibration characteristics of 
rotor blades was carried out. Coupled equations of motion for flapwise bend- 
ing and torsion were formulated for rotor blades with non-collinear elastic 
and mass axes. The finite element method was applied for a detailed repre- 
sentation of blade structural properties. Coupled structural mass and stiff- 
ness coefficients were evaluated. The range of validity of a set of coupled 
equations of motion linearized with respect to eccentricity between elastic 
and mass axes was. investigated. The sensitivity of blade vibration charac- 
teristics to torsion were evaluated by varying blade geometric properties, 
boundary conditions, and eccentricities between mass and elastic axes. 


INTRODUCTION 

A complete analysis of the dynamic response of helicopter rotor blades 
under actual flight conditions is difficult because of the three-dimensional 
nature of both blade deformations and unsteady flow around the blades. High- 
er operational forward speed and maneuverability requirements for modern hel- 
icopters have further increased the complexity of the problem. As indicated 
by theoretical and experimental studies, at higher tip-speed ratios the 
blade dynamic response characteristics become increasingly sensitive to tor- 
sional oscillation. 

The objective of the conducted investigation was to evaluate the effect 
of torsional oscillations on blade dynamic response characteristics and to 
analyze the sensitivity of torsional oscillations to blade aerodynamic, 
structural, and geometric properties. In the first part of the analysis, 
the coupled equations of torsional and bending motion of a rotating blade 
with non-collinear elastic and mass axes were developed in finite element 
form considering only structural parameters. The results from this formu- 
lation were compared with results from a system of differential equations 
of motion derived by Houbolt and Brooks for the same problem [1]. The im- 
portance of the second order terms in eccentricity between mass and elastic 
axes was established. The obtained overall results were also compared with 
results obtained by Isakson and Eisley [2], 
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a scalar 

generalized coupled mass and stiffness terms 
(see Eq. (A13)) 

a scalar 

a scalar 

2 

bending stiffness of a blade, N-m 
a scalar 

vertical displacement of a blade at the elastic 
center, m. 

column vector representing bending displacements 

vertical displacement of a blade at the center of 
mass, m. 

eigenvector of Eq. (B3) 

eigenfunction corresponding to the bending displacement 
component of p_. , m. 

2 

torsional inertia of the rotor blade, N-sec. 
identity matrix 

2 

torsional stiffness of a blade, N-m 
total stiffness matrix 
total stiffness matrix of an element 
stiffness matrix in bending 
coupled stiffness matrix 
stiffness matrix in torsion 

stiffness matrix in torsion for a blade element 
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M. 

—x 

Mb 

Me 

Mr 

Mr. 
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R 

s 

s. 


s, . 

ki 


t 
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U 

U E 

U„ 


length of blade element, m. 

2 2 

mass of the blade per unit length, N-sec. /m 

total mass matrix 

total mass matrix of an element 

mass matrix in bending 

coupled mass matrix 

/ 

mass matrix in torsion 

mass matrix in torsion for a blade element 
a scalar (see Eq. (A21)) 

column vector defining nodal forces in bending 
a scalar (see Eq. (A21)) 
rotor radius, m. 

2 

non-dimensional eccentricity, S = x^ ih/Iq 

non-dimensional eccentricity (see Eq. (A21)) 

non-dimensional eccentricity (see Eq. (A21)) 

generalized coupled bending damping in torsion, 

N-sec. /rad 

time variable 

centrifugal force along the longitudinal axis of a 
rotating blade, N 

total energy for coupled bending and torsion, N-m 
total energy in bending, N-m 
total energy in torsion, N-m 

distance on the longitudinal axis of cantilever blade 
from the fixed end, m. 

eccentricity between the elastic and mass center, m. 
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a scalar 

column vector representing bending and torsional 
displacements of a blade 

non-dimensionalized abscissa of a point on a blade element 

pitch change due to blade torsion, rad 

column vector representing torsional rotations 

torsional rotation of a blade at the i^ mode, rad 

square of the ratio of lowest torsional frequency to 
lowest bending frequency 

ratio of torsional rotations to bending displacements in 
the vibration mode corresponding to the lowest frequency, 
rad/m . 

column vector representing bending displacements 

coupled frequency of the blade, rad/sec. 

angular velocity of a rotating blade, rad/sec. 

natural frequency of a blade in bending, rad/sec. 

lowest coupled frequency corresponding to a predominantly 
bending mode, rad/sec. 

natural frequency of a blade in torsion, rad/sec. 

lowest coupled frequency corresponding to a predominantly 
torsional mode, rad/sec. 


Superscript 

t transpose 
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THE ANALYSIS OF STRUCTURAL VIBRATION CHARACTERISTICS 

OF ROTOR BLADES 

Method of Solution 

The energy method was applied in formulating the coupled equations of 
motion for helicopter rotor blades in flapwise bending and torsion. The 
formulation is similar to the one used by Garland [3] for analyzing non- 
rotating beams. The details of the mathematical analysis are given in Ap- 
pendix A. 


The coupled vibration characteristics of a rotating blade can be anal- 
yzed by first calculating the generalized coupled mass, damping, and stiff- 
ness coefficients in Eq. (A13) and then solving Eqs. (All) and (A12) . A 
coupled set of equations of motion, (Bl), was formulated as shown in Appen- 
dix B. A discretized set of mass and stiffness matrices was developed in 
finite element form. The derivation of each of the matrices in the set of 
matrix differential equations is given in Appendix C. 

A numerical method was applied for the solution in Eqs. (All) and 
(A12) . The following characteristics were established: 

* coupled structural vibration characteristics of rotor blades 
in bending and torsion 

• stability characteristics of the equations of Houbolt and Brooks. 


Variation of Coupled Vibration Characteristics with 
Eccentricity between Elastic and Mass Axes 

The variation of blade frequencies with increasing eccentricity between 
elastic and mass axes can be analysed from Eq. (A21) . This analysis is de- 
scribed in Appendix A. The variation of the frequencies of the sample blade 
for a two-degree-of-freedom system is illustrated in figure 1. As can be 
seen, when the non-dimensionalized eccentricity parameter S increases, the 
lower of the two frequencies corresponding to a predominantly bending mode 
decreases slightly, while the torsional frequency increases. 

A further consideration of figure 1 indicates the main features of the 
coupled structural behavior of a rotating blade, which can be summarized as 
follows : 

a) for values of the eccentricity Xq equal to zero (S = 0) un- 
coupled bending and torsional frequencies correspond to q = 1 
and q = A, respectively. As the coupling increases, the fre- 
quency corresponding to a predominantly bending mode decreases, 
while the torsional frequency increases. As can be seen from 
the left portion of the curve, the bending frequencies are 
bound between 
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2 2 
> a) > Q 


(1) 



b) the left portion of the curve can be approximated by a 
straight line as follows: 

q 1 - yS 


(2) 


The ratio of torsional rotations to bending displacements can be calculated 
from Eqs. (A18) and (A19) as 


2 

n 


/ VjA 

Vi 



2 2 


2 2 
S(fl -w ) 


m 



Observing also that 


2 2 
» W 

1 1 


and 


2 

n 


m 

I, 


2 

% Y 



S 


(3) 


(4) 


(5) 


the ratio of coupling between torsional and bending modes is 


n 


mu) r 


I 0 W T, 


( 6 ) 


Considering the left portion of the curve in figure 1, one can deduce from 
Eq. (6) that the coupling effect increases linearly with increasing eccen- 
tricity. 


Stability of the Set of Linearized Coupled 
Equations of Motion of Rotor Blades 

A set of equations of motion describing the vibrational behavior of a 
rotor blade has been developed by Houbolt and Brooks [1] in the following 
form: 
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Eqs. (7) and (.8) have been used by several investigators for calculating 
dynamic response of rotor blades [2,4]. These equations represent a linear- 
ized form of Eqs. (All) and (A12 ) with respect to x • only the first order 

terms of being considered. The second order terms E^ and FL correspond 

to the change in moment of inertia of the blade and change in the position 
of the axial force due to x Q , respectively. 

In this portion of the work the effect of linearization of Eqs. (7) 
and (8) was investigated in two steps: 

a) stability of the system of linearized Eqs. (7) and (8) 

b) comparative evaluation of approximations in Eqs. (All) and (A12). 

The sensitivity of frequencies to eccentricity between elastic and mass axes 
was calculated from Eqs. (7) and (8). Neglecting the higher order terms of 
x Q in Eqs. (A18) and (A19), these equations can be written in matrix form as 

follows : 



The dynamic system represented by Eq. (9) can be defined in terms of mass 
and stiffness matrices as 

[K - oj 2 M] 6 _ = [0] (10) 

The stability of the system represented by Eq. (10) can be determined 
by ascertaining whether the mass and stiffness matrices are positive defin- 
ite. For Eq. (1) one can show that 


2 

x m 

j — < 1 det [M] > 0 

0 


(ID 



det [K] > 0 


( 12 ) 
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Since 


C 13) 



> 1 , 


the rotor blade vibrations become unstable in the respective first torsion- 
al and bending modes, if 


2 

V 


> 1 . 


(14) 


The eigenvalues of Eq. (10) can be calculated from the following equation: 


det 






(A-q) 


(i-q) 


s 2 Cq-p) 2 


o 


(15) 


The non-dimensionalized eccentricity between mass and elastic axes 
related to the non-dimensionalized coupled frequency is shown in figure 2. 
Several conclusions about the coupled behavior can be derived from the 
curve shown in figure 2. From a comparison of Eq. (A21) and Eq. (15) it 
follows that 


= (l-q)CA-g) 
i (p-qHp-i) 


and 


q = (A-q) (i-q) 
2 (p-q) (p-q) 


From Eq. (16) one can derive the relation 


(16) 



(17) 


Thus, one can conclude that the solution of Eq. (15) is satisfactory for the 
left portion of the curve in figure 2. However, there is an important dif- 
ference for the right portion of the curve, since is not close to 

unity. 

The following conclusions can be deduced from the above results: 

a) the values of the lower coupled frequencies obtained from 
the Houbolt and Brooks equations would be satisfactory 

b) Eqs. (7) and (8) become unstable when there is a small 
increase in eccentricity between the mass and elastic axes. 

Such results are significant for the numerical integration of the equa- 
tions of motion. Instability of these equations leads to large torsional 
displacements and thus to the conclusion that the actual system is unstable. 
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Coupled. Vibration Analysis of a Sample Rotor Blade 


A sample rotor blade was chosen to illustrate the coupled vibration 
behavior. The finite element method, described in Appendices B and C, was 
applied for the numerical analysis. The structural properties of the sam- 
ple blade are listed in table 1. The uncoupled frequencies shown in table 2, 
were calculated using the numerical procedure summarized in Appendices B and 
C. The mode shapes for uncoupled vibrations are illustrated in figure 3. 

The coupled torsional behavior of the sample blade was then analyzed. 

Stability of the Houbolt and Brooks equations for the sample blade - 
The analysis described in Appendix B was applied for determining the stabil- 
ity of the Houbolt and Brooks equations; checking where the torsional modes 
become unstable. Considering only two-degrees-of-freedom for the sample 
blade, one can evaluate the critical eccentricity from 

2 A l°l 

Xq = , x Q = 0.089 (18) 

C 11 


For a four-degree-of -freedom system (two in bending and two in torsion), a 
quadratic expression was derived from the determinant of Eq. (Bll) for de- 
termining the critical eccentricity. The quadratic equation is 


A 1 A 2 D 1 D 2 " x e 


V C 22 D 1 * C 12 D 2’ * V C 11 D 2 * C 21 D 1> 


- X, 


C 11 C 22 ' C 21 C 12 


= 0 


(19) 


From Eq. (19) the value of the critical eccentricity is obtained as 

x 0 = 0.057 (20) 

The variation of coupled frequencies with increasing eccentricity was 
computed from the procedure described in Appendices B and C. The obtained 
numerical results for a two-degree-of-freedom system are shown in figures 4 
and 5 and listed in table 3. The variation of frequency with increasing ec- 
centricity agrees with the results predicted by Eqs. (16) and (17). Al- 
though the bending frequencies are quite accurate, the torsional frequencies 
can show considerable deviation. The instability of the torsional mode, as 
indicated in figure 5, agrees with the results from Eq. (18). For a multi- 
degree-of-freedom system the critical eccentricity that causes instability 
was found to be in the range of 



These results agree with the results obtained from the two-degree-of-freedom 
system in Eq. (19). Thus, for the sample blade the equations of Houbolt and 
Brooks may become unstable for moderate eccentricities. One has to be care- 
ful when using such equations as a numerical integration basis for the de- 
termination of the response of rotor blades. On the other hand, if desired, 
the second order terms in x^ can be included quite conveniently in the anal- 
ysis by using the finite element formulation. 

Effect of blade support conditions - The sample blade was analyzed 
with an assumed hinged support, which allows flapwise rotations. The cal- 
culated frequencies for the lowest modes and for various blade rotational 
speeds are presented in figure 6. The sample blade was also analyzed for 
a fixed-end boundary condition, in flapwise bending. The vibration fre- 
quencies for this case are presented in figure 7 and compared with the fre- 
quencies of a hinged blade as shown in figure 8. 

The coupled mode shapes of both hinged and fixed blades are shown in 
figures 10 and 11. As the end-fixity of the blade against bending rotation 
is increased, the bending stiffness increases; thus, the frequencies become 
correspondingly higher. The torsional component of the lowest coupled mode 
is also increased. 

Effect of eccentricity between mass and elastic centers - As shown in 
figure 7, the. lowest coupled frequency is affected by increasing the eccen- 
tricity between mass and elastic axes for the fixed blade; however, this 
frequency is less affected for the hinged blade. Frequencies corresponding 
to bending and torsional modes are compared for different eccentricities as 
shown in figure 12. Such results agree with the previously described re- 
sults on the variation of frequencies with increasing eccentricity. The 
torsional component of the lowest coupled mode is plotted for different ec- 
centricities in figure 12. The results also show the linear relationship 
between eccentricity and the degree of coupling observed in Eq. (6). 

Effect of rotational speed - The variation of frequency with increas- 
ing blade rotational speed for different boundary conditions and varying ec- 
centricities is shown in figures 6, 7, 8, and 12. As can be seen from these 
figures, the square of the coupled frequency term is linearly proportional 
to the square of the rotational velocity. 

The torsional rotation at the blade tip is plotted in figure 9 for 
various eccentricities at different rotational speeds. From the plotted 
curves one can observe that the torsional coupled motion increases linearly 
with eccentricity. For constant eccentricities this motion is also propor- 
tional to the square of the blade rotational speed. 

Effect of blade structural properties on coupled torsional behavior - 
It was shown that an approximate relationship, Eq. (20), existed for the 
coupled motion of rotor blades in bending and torsion. The validity of this 
relationship was examined by varying each of the parameters in the equation. 
The sensitivity of the torsional coupling parameter n with respect to each 
of the parameters was also checked for a multi-degree-of-freedom system. 
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When the distribution of the blade weight was uniformly increased by 
30%, it was found that the bending and torsional frequencies decreased; the 
first bending frequency by 1.5% and the first torsional frequency by 2.3%. 

In accordance with Eq. (A20) , the increase of the m/I ratio resulted in an 

0 

increased coupling of the response. As shown in figure 13 the torsional 
displacements increased by 30%. It appears then that the torsional coupling 
is quite sensitive to weight changes. 

Increasing the torsional inertia by 30% resulted in no change in the 
lowest bending frequency; however, the torsional frequency decreased as 
shown in figure 15. As shown in figure 14, there was also an increase in 
coupled torsional rotation in the second mode; however, there was no effect 
on the first mode. Such results can also be obtained from Eq. (A20). The 
increased torsional inertia and the resulting decrease in torsional frequen- 
cy has a balancing effect; producing no increase in coupling. The effect of 
decreasing torsional frequency is significant for the modes above the first 
bending mode and causes a net increase as shown in figure 14. 

The effect of an increase in the torsional rigidity on the frequencies 
and the mode shapes is presented in figures 16 and 17. Increasing the tor- 
sional rigidity by 30% led to a 25% reduction of the maximum torsional ro- 
tation in the first coupled mode. The lower coupled frequency was not af- 
fected, although the predominantly torsional frequencies increased. The 
plotted results correspond to Eq. (A20) ; showing an increase in torsional 
frequency and a decrease in the torsional component of the lowest coupled 
mode. 


Increasing the bending inertia by 30% resulted in no significant 
change in the bending and torsional frequencies or in the respective mode 
shapes. The results confirm again the validity of Eq. (A20) , since coupling 
is not affected by a variation of the bending inertia. 


Concluding Remarks 

The presented analysis extends the application of the equations of 
motion of Houbolt and Brooks for consideration of larger eccentricities be- 
tween rotor blade elastic and mass axes. The presented results indicate the 
significance of individual structural parameters in determining the sen- 
sitivity of coupled bending and torsional vibration characteristics. 

The results obtained for the sample blade show that for a rigorous 
analysis of vibration response characteristics, one should consider that: 

. the blade boundary conditions are important when determining 
the coupling ratio of torsional to bending displacements 

. the degree of coupling is linearly proportional to the square 
of the eccentricity between the elastic and mass axes 
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a more accurate representation of mass and stiffness properties 
is required for the coupled analysis of bending and torsion. 

When a detailed representation of the significant parameters is available, 
the blade sensitivity analysis can be used as a means of improving blade 
vibration and stability characteristics. 
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APPENDIX A 


DERIVATION OF THE EQUATIONS OF MOTION OF A ROTATING 
BLADE IN FLAPWISE BENDING AND TORSION 


The total energy of the blade, when coupling is not considered, con- 
sists of 


U B = 2 


"2 22 ’2 
(Elh - u>t mh + Th ) dx 

D 


(AI) 


for bending motion and 


U T= 2 


'R 


* 2 2 2 2 2 
(JG0 - «“ T 0 + a 10 ) dx 

10 0 


(A2) 


for torsional motion. 

For the coupled motion of the blade the displacements at the center of 
mass of the blade can be defined as 


h Q = h * x 0 e 

whereby Eq. (A3) is valid only for small values of 0. 

The potential energy for the coupled motion can be expressed as 
-R 

"2 22 »2 »2 22 22 
(Elh - (0 mht, + Th„ + JG0 - o>I Q 0 + n I ft 6 ) dx 

b b u o 


(A3) 




Substituting Eqs. (Al) , (A2) , and (A3) in Eq. (A4) , one can obtain: 
■R 




2 2 2 2 2 2 222 2 
(oig - a) ) mh + (w^ - u ) I Q 0 - a) m0 Xg - id m0x Q h 


2 *2 » » 

+ TXq 0 + 2Tx 0 h 
o y 


dx + Ug + U T 


(A5) 


For the bending and torsional displacements one can assume the fol- 
lowing functions: 
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00 


h = 


i=l 


a. h. 
1 1 


(A6) 


8 = 


j=l 


b.e . 
3 3 


(A7) 


where the displacement functions h^ and 0^ represent eigenvectors of the 
uncoupled Eqs. (AI) and (A2). One can also write 

'R 

T in o o i n 1 

(A8) 


„ _ 1 \ 2 
U B 2 / a i 


"2 22 '2 
Elh. - m„mhf + Th. 
i B l l 


dx 


and apply the principle of minimum energy to obtain 

3U t 


B 


3a. 


= 0 


(A9) 


For torsional motion one can similarly write 


3b. 

1 


= 0 


(A10) 


Substituting Eqs. (A6) and (A7) in Eq. (A5) , and considering the re- 
sults of Eqs. (A9) and (A10) , according to the principle of minimum energy 


a. (co 2 - ai ) A. - x b . 
l v B. J i 6i 

l 


2 2 
« B. . - o> C. - 
13 ij 


= 0 


(All) 


b j C<4 - » 2 ) D j - x e a i 


2 2 
n B. . - 0) C. . 

1 3 


- i/xl E. + n 2 xl H. = 0 (A12) 

9 3 0 J 


where 


A. = 

l 


-R 


mh. dx, 

l ’ 


-R 


B. . = d_ 
1J fi 2 


Te!h! dx, 
3 3 


14 . 





inh. 0 . 
i 1 


dx. 


E. = 
3 



D. = 
3 



Vj dx > 



(A13) 


C. . 
13 


9 


It is now possible to calculate the generalized coefficients 
D_. , Ey If ) and to investigate the effects of the coupling. 



B. . 

ij 


9 


An approximate investigation can be carried out when one considers the 
first bending mode and the first torsional mode of a rotating blade; i.e., 
for i = 1 and j = 1, Eqs. (All) and (A12) become 


“ 2) A i ‘ x e b i 


“Si - -Si 


= o 


V\ - » 2 ) Dj - x 9 aj 


“ 2b h - “ 2c n 


2 2 2 2 

“ x e E i + n x e H i = 0 


(A14) 

(A15) 


The lowest uncoupled modes in bending and torsion, respectively, of a 
rotating blade of uniform cross section can be approximated as 


h l = e l x 

(A16) 

e l - f l x 

(A17) 


For a rotor blade of uniform cross section and uniform eccentricity 
along the blade, Eqs. (A14) and (A15) become 


2 2 2 2 
(u) - u) ) ma.e + (SI - co ) mx Q b 1 f. l = 0 

11 D 1 1 


(A18) 


2 2 

(S2 - m ) mx 0 a 1 e 1 + 


(<4 - a) 2 ) I 0 + mx 2 (SI 2 - co 2 ) J> b 2 f 2 = 0 (A19) 


From Eqs. (A18) and (A19) we can now write 
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The expression (A20) can be represented as 

s = (q-1) (q-A) 
1 (p-q) Cp-1) 


(A20) 


(A21) 
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APPENDIX B 


FINITE ELEMENT FORMULATION OF THE LINEARIZED EQUATIONS OF 

HOUBOLT AND BROOKS 


A general solution to the equations of Houbolt and Brooks can be ob- 
tained by formulating the equations in matrix form and then applying the 
finite element theory. These equations can be written in general form as 
follows : 


/ 

- 




x M 


H3 

0— C 

2 



- u 




wt 



x e“c 

— T 


— 

— 1 





> 


r 

K 

x 



h 

— B 

0-^ 







> 


_ X e^ 

*T_ 

/ 


0_ 


(Bl) 


The formulation of the matrices in Eq. (Bl) is given in Appendix C. Total 
energy of the blade can be expressed in matrix form as follows: 


.. 1 2 
U = - j a) 


^Mgh + XgQV^h + x 0 h t M c £ + 


1 

+ 2 


t 

— — B^- 


V Q— C- 


h K„h + xj Ch + x^e + 


(B2) 


For the displacement vectors h and 0_ one can assume expressions of the 
following form: 


2 

V 

i 


♦JMk - 

[0] 

h = a.h. 

— l—i 

2 

“t. 

1 

^r-j 

II 

®r 

+ 

[0] 

0 = b . 0 . 
- J-J 


Eq. (B2) can then be written as 


U = - i- w 2 a 2 A. + 2x a. b .C. . + b 2 D. 

2 [ i i 0i)i] 1 1 J 

+ ^ a?h*]Ui. + 2x a. b . B. . + b 2 0*K_,8. 
2 i— i— T — i 9i]i] J — J~T" _ 1 


where 


A. = hu M„h . , 
l — l— B— l 

B. . = e*K_h. 
ij -i-c-i 

c. . = e^MJi. , 

il -j-C-i’ 

o 

II 

Jf 0 


(B3) 

(B4) 


(B5) 

(B6) 

(B7) 


17 



An application of the principle of minimum potential energy on Eq. (B5) to- 
gether with Eqs. (B3) and (B4) leads now to the following system of equa- 
tions : 


"2 2 


2 

to_ - to 

a. A. + x.b. 

B. . - to C. . 

B. 

l 

ii 6 j 

il ij 


(B8) 


'2 2 

b.D. + x„a. 

2 

to™ - to 

B. . - to C. . 

T . 

L i J 

11 9i 

il il 


(B9) 


Eqs. (B8) and (B9) can be written in matrix form as follows: 


(tOg i -to 2 )A 1 

I V B ir“ c n ) • • 


1 

h-» 

1 


1 

o 

1 

(to^-to 2 )^ 

' x e (B 2 r“ 2c 2 i ) 


a 2 


0 

• 

• 


• 


• 


1 

• • 


a 

m 

b i 


• 

x e CB n-“ 2c ii 3 

J- 

• | t“T 1 -“ 2)D 1 


= 

0 

X 0 ( - B 12“ a) C 12' ) 

1 (4 “ w2 ) D 2 

| 


b 2 


0 

• 

1 


• 


• 

• 

l 


b 

n 


• 

_ 

1 


„ _ 


_ _ 


(BIO) 


The terms <o n , to™ , A., B. . , C. and D. in Eq. (BIO) can be evaluated for 
B. T. 1 ij il l n ' 

l j J 

an uncoupled system of equations. Eq. (BIO) can then be solved to obtain 
coupled frequencies. 


The stability of Eqs. (B8) and (B9) can be tested by determining 
whether or not the mass matrix is positive definite. In matrix notation 
this can be shown as 
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APPENDIX C 


FINITE ELEMENT FORMULATION OF MASS AND STIFFNESS MATRICES 
FOR THE COUPLED EQUATIONS OF MOTION 


The finite element method was used for a discretized representation of 
the mass and stiffness properties of the rotor blade. The coupled equations 
of motion (Eqs. (7) and (8)) were written as a system of matrix differential 
equations. To analyze the torsional vibrations, a constant strain element 
was used which requires a linear variation of displacement over the element, 
while for bending vibrations a quadratic displacement function was employed. 


Formulation of Mass and Stiffness Matrices for Torsional Vibrations 


A consistent mass matrix was used for the discretized representation of 
the inertia term in the equation of motion (Al). The derivation of the mass 
matrix for torsional vibrations will now be summarized. For the blade ele- 
ment in figure 18, with nodal torsional rotations 0^ and 0 , the rotation 
at any point can be expressed as ' 


6 = [1-C C] 


(Cl) 


The torsional acceleration at any point can then be written as 



C] 

1 

a 

K> 

CD 

— -[1C 
9 L 1 

dt 

dt 2 



d 2 e 



7 1 


(C2) 


The equivalent inertia forces at the nodes of a blade finite element can 
be defined as a vector T using the principle of virtual work as follows: 


0 t T = 0* 



1-C 


[i-c e] dx 0 


CC3) 


or 



2 


1 


T = 


T e £ 



The elemental torsional mass matrix can be written as 



: e * 

6 


2 

1 


1 

2 


(C4) 


(C5) 


The elemental stiffness matrix can be analyzed in two parts. The first 
part is the elastic stiffness matrix, while the second part is the matrix 
representing the effect of the centrifugal forces along the blade axis. For 
a linear variation of torsional displacements over a finite element, using 
an approach similar to the derivation of the mass matrix, the effect of the 
centrifugal forces can be represented by a geometric stiffness matrix. The 
torsional stiffness nratrix for a rotating blade element can then be written 
as 



+ Q 


2 *0 A 


1 

2 


(C6) 


Formulation of the Mass and Stiffness Matrices for Bending Vibrations 


After defining the nodal force vector £ and the nodal displacement vec- 
tor p_, as shown in figure 18, the vertical displacement, h, at any point 
can be expressed in terms of the nodal displacements as follows: 


(1 - 3S 2 + 2C 3 ) *(c - 2c 2 + c 3 ) (3C 2 - 2C 3 ) * C 3 ) 


(C7) 



A consistent mass matrix can also be derived for the bending vibrations of 
the finite element using the approach applied for torsional vibrations. 
Neglecting the effect of the shear deformations, the mass matrix in bending 
can be written as 


21 



r 


156 





, 2 

Symmetric 

M„ = 

m & 

222 

42, 

— B. 

l 

420 

54 

132, 

156 



- 132 

-32, 2 

-222, 4 2, 2 


CCS) 


In the case of a rotating blade the bending stiffness can also be analyzed 
in two parts as elastic stiffness and geometric stiffness due to the axial 
forces in the blade. The representation of the elastic stiffness of a blade 
element is well known [5], The effect of rotation of the blade can be eval- 
uated by making use of a geometric stiffness matrix* which has already been 
used for solving stability problems [5], 


From Eq. (Al), the following relationship can be written for the finite 
element to represent the virtual work by axial forces: 



h Th dx 


o 


t 


(C9) 


Substituting Eq. (C7) into Eq. (C9), the bending stiffness of the blade ele- 
ment due to axial forces can be calculated. The total stiffness matrix in 
bending can then be written as 



12EI 





6 




2 3 





5 




6EI 

2 2 

1 

so tn 

1 1—1 

Symmetric 


2, 

10 

22 2 

15 

Symmetric 





T 

+ — 





-12EI 

-6EI 

12EI 


2 

-6 

-2 

6 


2 3 

2 2 

2 3 



5 

10 

5 


6EI 

2EI 

-6EI 

4EI 


2 

-2 2 

-2 

2 2 2 

2 2 

2, 

2 2 

2, 


10 

30 

10 

15 


(CIO) 


The displacement functions 
Eq. (C7) , and for torsion, Eq. 
the coupling terms in Eqs. (7) 
placement functions. 


were defined fora blade element in bending, 
(Cl). Using the principle of virtual work, 
and (8) can be derived for the defined dis- 
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Derivation of the Coupled Mass Term (mx 6) 

The virtual work due to the coupled mass term can be expressed as 

t P 

p_ P = | h m x Q 0 dx (Cll) 

Jo 

.Substituting Eqs. (C2) and (C7) into Eq. (Cll) one can obtain Eq. (C12). 

21 9 

3£ 2£ 

9 21 

-2£ -3£ 

' 2 » 

Derivation of the Coupled Stiffness Term (ft mr x ) 

0 

The virtual work due to this term can be formulated as 




2 1 

ft mr x. 0h dx 
0 


Assuming that the length of the finite element is small in comparison to r, 
Eqs. (Cl), (C7) , and (C13) yield the following expression: 



Determination of the Coupled Mass and Stiffness Matrices 

Mass and stiffness matrices in bending and torsion can be combined to 
obtain a single mass-stiffness representation of the finite element for 
describing the coupled behavior. Combining the mass and stiffness matrices 
from Eqs. (C5), (C6), (C8) , (CIO), (C12), and (C14), and defining the 
force and displacement vectors for the element as 
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V— 1 

P = 

" P l‘ 

P 2 


P 2 

T 1 


CD 

P 3 


P 3 

P 4 


P 4 

T 2 


9 2 


(C15) 


the coupled mass and stiffness matrices for the element can be written as 
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4£ 2 
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Symmetric 
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63Xg 

14£Xg 
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TABLE 1. CH-34 Blade Physical Properties: 
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TABLE 2. Uncoupled Frequencies of CH-34 Blade in Bending and Torsion: 


Mode No. 

Bending (rad/sec.) 

Torsion (rad/s 

1 

25.2 

173.6 

2 

72.4 

532.1 

3 

136.2 

934.9 

4 

210.8 

1350.0 


TABLE 3. 


The Frequencies Corresponding to the Coupled Modes, Consisting 
Mainly of First Mode in Bending and First Mode in Torsion. 


x. (m) % (rad 2 /sec. 2 ) u> 2 (rad 2 /sec. 2 ) 

0 B 1 


0.0 

635.0 

30,137 

0.020 

632.8 

31,959 

.041 

626.2 

38,499 

.050 
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47,284 

.060 

616.3 

59 , 349 

.069 

611.0 

83,710 

.077 

604.9 

146,223 

.083 

600.5 

289,632 

.104 

581.9 

-71,923 

.124 

559.6 

-31,572 

.145 

539.0 

... _ -18,779 

.187 

473.7 

-9,520 
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TABLE 4a. Generalized Structural Mass M £.), Xg = .0254 m for Uncoupled Modes 

of CH-34 Blade 
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TABLE 5. Uncoupled Generalized Structural Mass (ol), Critical Damping (3^) > and Stiffness (y^) 
of CH-34 Blade 
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to a Mainly Bending Mode with Increasing Eccentricity. 



Frequency Ratio - q 

Fig. 5 Variation of the Lowest Frequency Corresponding to a 
Mainly Torsional Mode with Increasing Eccentricity. 
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1st Mode 



Fig. 6 Variation of Lowest Bending Frequencies for a Hinged Blade 
fCH-34) with Blade Rotational Velocity, (x. = 0.0 and 3.28) 
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Fig. 7 



Variation of Lowest Bending Frequencies, for a Cantilever 
Blade (CH-34) with Blade Rotational Velocity. 
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Fig. 8 Variation of Lowest Bending Frequencies for a Cantilever 
and Hinged Blade (CH-34) with Blade Rotational Velocity. 
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Fig. 9 Variation of the Torsional Rotations at the Blade Tip 
with Eccentricity. 
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Fig. 12 Interaction Between Mainly Torsional Modes and Higher Bending 
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Fig. 16 Variation of First Coupled Mode with 30-6 Increase in 
Torsional Rigidity (CH-34). 
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Fig. 17 Variation of Torsional Frequencies with 30% Increase in 
Torsional Rigidity (CH-34) 





Force Displacement 

Fig. 18 Nodal Forces and Displacements for a Blade Element in 
Bending and Torsion. 
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